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Abstract 


Using lattice QCD on a 16^ x 24 lattice dX j3 = 6.0, we examine the 
elastic limit of charge overlap functions in the quenched approximation for the 
pion and rho meson; results are compared to previous direct current insertion 
calculations. A good signal is seen for the pion, but the electric and magnetic 
rho meson results are considerably noisier. We find that the pion and rho 
results are characterized by a monopole mass to rho mass ratio of 0.97(8) 
and 0.73(10), respectively. Assuming the functional form of the electric and 
magnetic form factors are the same, we also find a rho meson g-factor of 
g = 2.25(34), consistent with the nonrelativistic quark model. 
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I. INTRODUCTION 


The spacelike pion form factor has been the subject of a number of completed lattice 
studies However, it has not yet been studied systematically by the method of Ref. 

It would be helpful to compare the systematic and statistical characteristics of current 
overlap to previous results as well as to learn about possible limitations. Since the same set 
of quark propagators can also be used to construct current overlap functions for higher spin, 
we investigate the rho meson as well. As was the case with the decouplet baryons studied 
in Ref. |^, the rho meson is not stable under the strong interactions. However, at quark 
masses available in lattice calculations, it’s decay is kinematically forbidden, so in this sense 
our quenched results are physical. Another motivating factor in studying the rho is to have 
a more “typical” hadron to compare our pion results to, especially in terms of the relative 
statistical signals. Our results can also be used to compare to hadronic models. 

Current overlap techniques are versatile in that the same set of quark propagators can be 
used to study both elastic and inelastic processes. In the second part of this series, referred 
to as Part H, the inelastic part of the correlation functions calculated here will be used 
to partially address the question of the polarizability of the charged pion. Thus, another 
purpose of the present study is to set the elastic “baseline” necessary to extract additional 
nonelastic properties of hadrons. In addition to polarizability, structure functions may also 
be extracted using these techniques [@,^. 

In summary, the purpose of this study and others like it is both as a testing ground for 
lattice techniques and as a preliminary contact between experiment and fundamental theory. 
As we explore the numerics, we will be more interested in trends in our lattice data rather 
than the final numbers. Nevertheless, we do not neglect comparison wherever possible to 
previous lattice and experimental results. We will partially explore the systematics associ¬ 
ated with longer time correlations, larger volume and finite lattice spacing effects. However, 
further work treating systematic effects such as scaling and quenching will need to be done 
before our numerical results can be accepted as “physical”. 
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We will start with a review of the formulas used in this study and then proceed directly 
to the results. We hnish with a summary and some brief comments about the possible 
physics underlying our results. 


II. CORRELATION FUNCTIONS 


A. Formulas 


The formulas necessary for measuring the form factors of the pion and rho meson by 
the current overlap technique have been developed in Ref. ||^. The relevant results are the 
following. One can dehne the u, d flavor charge and current density overlap matrix elements 
(continuum, imaginary time formalism) for these two particles as : 

^ ^(^+(o)|T[-/(r,()p"(0)]|ir+(O)), (1) 

= ^(p+(e;0)|r|V(f,i)p“(0)||/K;0)), (2) 

MC.Pt) = 5^(pyf;O)|r|V(f,i)J*"(0)||p+K;O)). ( 3 ) 

t) represents the charge density operator ( 0 ) represents the vector 

particle state of spin polarization .^ (.^ = ±, 0 refers to spin component along the z-axis) and 
zero momentum. We dehne the Fourier transforms: 



J dVe“-7p(f,t), 

(4) 

III 

(M 

'd3re“-7p(^;r,t), 

(5) 

Kk{^]q,t) = J 


(6) 


In the large Euclidean time limit we can design form factor measurements from the ampli¬ 
tudes (“a” is the lattice spacing) 




AEqUl. 




AEqrUp 


(7) 

( 8 ) 
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(9) 

-Gj(±;4,)e-l®--"-)‘, 

(10) 

) 

(11) 

can show that Gp{±]ql) 

= G'p(0;g2y). In 


Eq.fpTD we may cyclically change the directions of the current, polarization and momentum; 
anticyclic changes produce a minus sign. We can then extract the charge (Gc), quadrupole 


{Gq) and magnetic {Gm) form factors of the rho meson from (the common magnitude of 
spatial momentum is q) 


Gc — Gp(±; ql){- + -if), 




q^ 


2g,(+;a 

" G,(±;,|)(l + ,,)' 


where 


>) = (2 


G?(±; ql,) 

ayqi) 



( 12 ) 

(13) 

(14) 


(15) 


B. Simulation Details 

Our quenched lattices were constructed by the algorithm in Ref. |§], thermalized by 
11000 sweeps, and separated by 1000 sweeps. (Ten of these lattices overlap with those used 
in Ref. 0.) Our results are obtained with Wilson fermions on twenty 16^ x 24 lattices at 
f = 6.0. We used the exactly conserved lattice charge and current densities in forming the 
overlaps. For the quarks, we used periodic boundary conditions in space and “hxed” bound¬ 
ary conditions in time (lattice time boundary gauge held time links are unused). When 
a particular polarization and momentum state is called for in the above, we measure this 
in all possible ways in a given conhguration in order to reinforce the signal. For example. 
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the Gp{^] ql) amplitude can be measured by averaging the results for ^ = ±. An addi¬ 
tional two measurements are afforded by Gp{0]q‘^y). When a correlation function for the 
second momentum is constructed, we combine signals from q = (±|,±|,0), (±|,0,±|) 
and (0, ±|, ±|). 

In forming these amplitudes on the lattice, we neglect the charge and current self¬ 
contraction loops, which do not vanish in this context. In terms of quark lines, these are 
the disconnected (sea quark) diagrams, which are extremely difficult to simulate^. How¬ 
ever, it is unlikely that such local objects will signihcantly affect the elastic (large time 
separation) limit studied here. The connected amplitude and an example of a disconnected 
contribution are shown in Figs. 1(a) and (b), respectively. In collecting data from our prop¬ 
agators, we do so in a manner which is as symmetrical as possible in time, given that the 
conserved lattice charge density is nonlocal. Thus, we use the symmetrized charge den¬ 
sity, |[p(r, t) -|- p{f,t — 1)], naturally associated with integer time locations, in forming the 
Hp{^] q ^) correlations. In addition, we hnd it is important to locate the two currents as 
symmetrically as possible between the particle interpolation helds, which are hxed in time. 
When a symmetrical time array of currents is not possible, we average over the two possi¬ 
bilities on each conhguration. The time behavior of our signals is then much smoother. 

Our lattice interpolation helds for the pion and rho are standard, given by and 

where e(±) = e(0) = i (rho meson rest frame). However, our inter¬ 

polation helds are smeared over the entire lattice spatial volume (using the lattice Coulomb 
gauge), projecting onto zero momentum. Form factors are calculated at three values of the 
hopping parameter {k = 0.148, 0.152 and 0.154) and parameters are extrapolated linearly 
to the chiral limit. 

The lattice analog of the coordinate space functions P{f,t), P{^]f,t) and 
which we shall denote as and can be shown to be real and 

purely imaginary, respectively, in the conhguration average. That is, using the identity 
developed in Ref. []^ and given that the gauge helds U and U* appear with equal weight 
in the ensemble average, one may show that 
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PK;r.t; {t/}) = T>K; f, i; {(/•})*, 


( 16 ) 

( 17 ) 


and 


r, t- {U}) = f, t- {[/*})*. (18) 

These identities allow us to neglect terms which are purely noise in the discrete Fourier trans¬ 
forms of V{f,t]{U}), V{i]r,t]{U}) and (P(r,f), P(^;r,f) and Jk{i]r,t) 

are clearly bad notation for these four-point functions which obviously depend on more than 
just the relative spatial and time separations. However, Eqs. ([T^) -(|l8D convey the basic idea.) 
Also, by combining Fourier transform measurements involving both q and —g as we do, one 
is explicitly symmetrizing or antisymmetrizing these, so there is no need to assume evenness 
or oddness in r. We can actually show that the identities. 


V{r,t) = 

(19) 

■,yt) = v{y,-yt)y 

(20) 


hold exactly, configuration by configuration, for the lattice charge overlap functions. We 
can find no corresponding exact statement for J^k{^]T,t), however. 

One important aspect of forming the overlap functions is to have an efficient algorithm 
for sewing together the quark lines to form the charge and current densities at all relative 
spatial separations. Our analysis codes use the following trick to save computer time. The 
inner loop has a sum that looks like: 

/(O = II/2((? + r)/l(f). (21) 

X 

(x, r are spatial locations in a cubic 3D space.) This loop, which scales like iV^, where W 
is the number of space points in the lattice, and must be repeated at all chosen relative 
time separations, is not easily vectorizable for spatially periodic boundary conditions. This 
equation was replaced by a version using Fourier transforms: 
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( 22 ) 


where 


/(,-0=i:e-«-’'Q2(-,-)Ql(j), 

<? 


=i:e-‘’"*71(Si). 

xi 

Q2(-g) = 


(23) 

(24) 


This version is easily vectorized (i.e., the Q2{—q) times Ql{q) above) and can ntilize fast 
Fonrier transform rontines^. This techniqne allows us to extract, with only a small amount 
of computer time, the full time extent of the current overlap correlations, even up to t = 22, 
where the currents begin to overlap with the interpolation helds. 

In making our choice of hts to masses and correlation functions, we examine the chi- 
squared per degree of freedom, using correlated hts. The covariance matrix |^, Cij, is 
estimated from the single-elimination jackknife 


N — 1 ^ 

C,1 = E(V(n) - A'i)(X,(n) - X,), 

n=l 


N 


(25) 


where N is the number of conhgurations, i and j label different time slices, and the Xi{n) 
represent jacknifed propagator data, being the average to remove the bias. The is 
then given by 


x"=i:cy(r. - /i)(y - *>. 


(26) 




where the Yi are the average experimental values of the propagators and the /, are the 
functional form values for these time slices. However, it has been noticed by many authors 
that the covariance matrix overestimates correlations on propagators for small numbers of 
conhgurations [Q. Therefore our procedure on hts is the following. We choose our time 
intervals based upon obtaining acceptable values of the Xd- However, following the suggestion 
in Ref. [Q, the actual hts themselves are uncorrelated, which simply means neglecting the 
oh-diagonal components of the Cij. We feel it is necessary, and will see later it can be 
important, to take correlations into account in extrapolating our results across k values. 
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Therefore, the full forms of Eqs.(^5D, (^6[) are used in our chiral extrapolations {i,j then 
label K values). 


A third order single elimination jackknife was used for error analysis; the hrst order 
dehnes error bars on the time correlation functions, the second dehnes error bars on the 
time correlation fits and the third is necessary for the chiral extrapolation of the results. 
Our masses were measured on twenty gauge conhgurations using single exponential fits and 
are listed in Table I along with the chi-squared per degree of freedom, of the £t. The 
source position for these propagators is the hrst time slice of the lattice (t = 0). The masses 
are consistent with those measured with twelve conhgurations in Ref. 0 . However, in order 
to hnd a time interval on which the correlated Xd acceptable for our pi and rho masses 
simultaneously, we had to increase the time interval considerably beyond the hts in this 
reference, resulting in slightly increased error bars. 

Two dehnitions for the Wilson quark mass, which agree for k ^ Kc, have been used in 
previous studies of form factors. These dehnitions are 


ma = 

2 K 



(27) 


and 


ma = Ini—- — 3), 

K 


(28) 


where Kc is the critical k value for which the pion mass vanishes. The hrst dehnition is 
motivated by leading order chiral perturbation theory; the second dehnition comes from 
tadpole improving |]^ the quark pole mass in free held theory. Fitting the square of the 
pion mass in Table I with Eq. (|^) gives Kc = 0.1566(2). (Ref. |T^ gives Kc = 0.1570(1)±.0002 
from f = 9 to 14 point-source hts, whereas the hts here are f = 15 to 18 smeared-source. 
See the next section for more comments on our time ht choice.) Fitting the same data to 
Eq.(p8D gives Kc = 0.1564(2). We prefer Eg. (p8|) and the latter Kc value in this work to 
partly correct for the fact that our quark masses may not be in the leading order chiral limit 


range. 




As was done in Ref. [Q, we shall concentrate here on dimensionless quantities in extrap¬ 
olating to the chiral limit. This is because a ratio of similar physical quantities is often less 
subject to systematic errors. Our philosophy in comparing our results to experiment is to 
assume the simpliest possible functional form and to extrapolate the £t parameters rather 
than individual form factor values in order to make contact with phenomenology. This is 
crucial in the case of the pion, where the four momentum transfer in the chiral limit vanishes. 

III. RESULTS 

A. Preliminaries 

First, let us make a point about our correlation functions, which are four-point functions. 
In order to measure current overlap functions, it is necessary to form an amplitude which 
looks generically like Fig. 1(a). Eqs. 0-0 tell us the form factors are identihed in the large 
time separation limit, t, of the charge or current densities. However, large time separation 
between the currents produces small time separation between the hxed interpolation helds 
and the currents. It is only as we increase the time separations between all four time locations 
of our four-point function that the lattice amplitudes project with increasing accuracy on 
the ground state. This is quite different from generic two-point functions which, outside of 
time boundary effects, are guaranteed to have a better ground state overlap for larger time 
separations. This has important implications for our time hts, as we will discuss below. 

We next discuss two tests done to help motivate the choices made in this study regarding 
the use of smeared helds and for the time position of the interpolation helds. 

One test is given in Fig.2(a), which shows a graph of the Fourier transform of the equal¬ 
time charge overlap function, Q(g^, t), as both charge densities are moved in unison between 
pion source and sink aX k = 0.154. This measurement is associated with half-integer time 
steps since the conserved lattice current is non-local in time. We are using smeared-to- 
smeared quark helds and are examining the two lowest momenta on a linear scale using 
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twenty configurations. The data should be flat if the lattice is long enough in time. The 
satisfactory results seen in Fig.2(a) are to be contrasted with those in Fig.2(b), which repeat 
the same measurement (using a subset of ten conhgurations), but using smeared-to-point^ 
correlation functions. 

Having adopted smeared helds, in Fig.3 we show measurements of the Fourier transform 
of the pion overlap distribution functions, Q((f^, f), on a logio scale at the two lowest spatial 
momenta transfer for AT = 17 (boxes) and AT = 23 (diamonds), again at k = 0.154. AT 
refers to the time separation of the pion source and sink helds; time separations between 
the charge densities range from 0 to 10. (The AT = 17 time correlation data is from 
Ref. 0 , where only charge overlap time separations up to 10 were considered because of the 
expensiveness of constructing the overlap function. See Section 11(B) above for comments 
on the improved method used here.) Note that the AT = 17 measurements use twelve 
conhgurations and the AT = 23 use twenty, ten of which overlap. There are no signihcant 
diherences within statistical errors for the lowest momentum (upper set of points), but 
there appears to be a diherence larger than the error bars on the second momentum results. 
To examine the situation from another point of view, consider Figs.4(a) and (b). Here 
we have plotted the local {E — m) determined from the exponential falloh of the pion 
overlap function for the three k values in this study in (a) for \qmin\ = f and (b) for 
1^ = \/2\qrain\- This measurement is also naturally associated with half-integer time values. 
Although the diherent k values approach the asymptotic limit at different rates, all of the 
hrst momentum results are quite consistent with the expected exponential falloff (determined 
from continuum dispersion and the measured pion mass) by a time separation of 14 steps. 
The second momentum results in Fig.4(b) are more problematical. The k = 0.154 local 
{E — m) results are seen to eventually approach the correct asymptotic falloff (slower than 
for the hrst momentum), albeit with large error bars. However, the smallest k value at this 
larger momentum is apparently displaying a violation of continuum dispersion. (The dotted 
lines in Figs.4(a) and (b) give the spin 0 lattice result; see Eq.(^) below.) Because of this 
apparent violation we prefer in this study to exclude higher momentum measurements for 
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both the pion and rho meson. (We will estimate the remaining effect of finite lattice spacing 
on our first momentum extrapolated quantities at the end of the next section.) An analysis 
of the lowest momentum rho meson local electric {E — m) shows a similar pattern to the 
pion with, however, larger errors. In addition, we will see later that the rho meson magnetic 
(E — m) becomes consistent with continuum dispersion much earlier in time than the electric 
ones. Thus, our lattice appears to have enough time steps to filter out asymptotic limit of 
the first momentum in the overlap correlation functions and continuum dispersion is reliable 
within errors. 


B. Numerical Results 

After these preliminaries, we now begin with a survey of the correlation functions mea¬ 
sured in this study, given by the quantities in Eqs.(|7|)-(pd]). The different n values are 
combined and shown in Figs.5-8. We will take some pains to examine all the correlation 
functions since this is the first complete study using the current overlap technique. 

First, let us explain our fitting procedure. As pointed out above, large time separations 
of the currents move them toward the interpolation fields, fixed at the time ends of the 
lattice. There is no d priori reason that the large time separation amplitudes are to be 
favored; non ground state contributions are expected both for small as well as for large time 
separations. We will use the calculated Xd of fo determine allowed fit time intervals 

for the correlation functions. In our fits, we wish to extend the time plateaus sufficiently 
that we are testing the fit in a nontrivial way, but not so extended that the becomes 
unmanageable or the statistical signal has decayed. Our choice of fitting 3 time values 
on overlap functions and 4 on mass correlation functions is thus a compromise between 
systematics and statistics. (Both types of fits involve 2 independent degrees of freedom.) 
Our overlap fits will assume the continuum dispersion relation, E'^ = q'^ + m?. We attempt 
to fit all K values for given type of correlation function on the same time interval in order to 
increase possible correlations across k. We find this be be a significant effect for the chiral 
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extrapolation of the pion data. (See the discussion under Fig.9 below.) 

Fig.5 shows the correlation functions, versus time separation, t, for the pion. 

(The next 4 hgures are logio plots.) The lines drawn are the best hts over a common f = 14 
to 16 time separation interval. A single exponential form is evident even for large time 
separations. The form factors, F, from these fits using Eq.(^ are listed in Table II along 
with the implied F / F^^ ratio {F^^ is the vector dominance result predicted by the masses 
in Table I) and the x% Notice that all of the F/F^^ ratios are low compared to the vector 
dominance result of 1. 

Fig.6 shows the time separation correlation function for the rho meson corresponding to 
adding the signals for Q(±;g^,f) and Q{0-,qly,t) (see Eqs.(8) and (9).) The error bars are 
signihcantly larger than for the pion. The k, = 0.148 and 0.152 results are £t across t = 14 
to 16 (as for the pion) with acceptable Xd values; see the Xd values in Table IV under the 
line for Gc/G^^. However, a similar fit to the k = 0.154 data gives a rather large value 
for Xd (1-99) caused by poor exponential behavior near the time separation edge. We have 
therefore moved this £t inward until an acceptable £t is achieved. This occurs after shifting 
it a single time step (t = 13 to 15 £t); we now have Xd = 0.43. We note that the local 
{E — m) values for k = 0.154 are low, but consistent within large error bars with continuum 
dispersion across these time slices; the other two k values are also consistent in this sense. 

Fig.7 shows the magnetic correlation function iFy(+;ga.,f) as a function of time separa¬ 
tion. The time behavior of these functions, which come from charge-current overlap (see 
Eqs.(j|) and (|^)), is quite different from the electric (charge-charge overlap) functions just 
examined. We are able to fit these functions much earlier in time to the required exponen¬ 
tial behavior with acceptable values than in Figs.5 or 6, resulting in improved error bars. 
The situation is similar to the magnetic correlation functions for the proton, examined in 
Ref. [^, which also assumed their expected exponential behavior more quickly than their 
electric counterparts. The values of these hts are listed in Table IV (under the line for 
Gm)- Again, the local {E — m) values across the chosen time slices are consistent within 
error bars with continuum dispersion. 
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Fig.8 shows the rho correlation function Q{±-, ql y,t) as a function of time separa¬ 
tion. There seems to be a nonexponential systematic effect for larger time separations 
for K = 0.152 and k = 0.154. In addition, the local {E — m) plots for these k values are 
never consistent within error bars with the expected exponential falloff. We feel that the 
appropriate conclusion from these difficulties is that the quadrupole form factor can not be 
reliably extracted with our data sample. We note that if the quadrupole form factor is zero, 
the results in Figs.6 and 8 should be the same within statistical errors for each k value; 
there is indeed strong overlap of the error bars. Therefore, in what is to follow, we assume 
that Gq{q^) = 0, consistent with the data^. In Eqs.(|^-(|T^ this makes rj = 1. The values 
for the electric, Gc, and magnetic, form factors in Table IV were extracted with this 
assumption. 

Now that we have examined the correlation functions, in Fig.9 we begin to show results 
of extrapolations of the pion correlation functions. The form factors from Table II are 
characterized by the corresponding value of the monopole mass, mM, implied by < 0) 


= (1 






(29) 


The ttim values arising from the form factors are divided by the measured rrip and displayed 
in Table III and Fig.9. Also shown are pion results from Ref. [§, which used 28 gauge held 
conhgurations of size 24^ x 12^ including an extra quark propagator with reversed momentum 
to decrease error bars. (Values and error bars for the pion and rho meson are extracted from 
the charge radius results in Table 1 of Ref. 0 using uim = VQ/r.) The result of our analysis 
is linearly extrapolated to the chiral limit; the Xd 0.53. Our hnal result for the mM/mp 
ratio is 0.97(8). Note that the error bars on rup are not included in the our results; this is 
not the conservative assumption. However, one reason for this is so our error bars can be 
compared with those in Ref. 0. In addition, uim and nip are likely to be strongly correlated 
in both studies, so niM/mp is probably better determined than uncorrelated addition of the 
error bars would indicate. Also note that the physical of Ref. [Q, set by the nucleon mass. 


is |g^| ~ 0.16 GeV^. If we set our scale in the same fashion using the results of Ref. [jT6 
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our physical momentum is \q^\ ~ 0.47 GeV^, approximately three times larger. Ref. 0 
uses (3 = 5.9 as opposed to 6.0 here which means that the quark mass times lattice scale 
values, ma, from Ref. 0 must be multiplied by the ratio of scales so that we can compare 
results at the same physical quark mass value. (Ref. § uses Eq.(^ for ma.) The lattice 
scale in Ref. is reasonably consistent within large errors with Ref. using asymptotic 
scaling. The nucleon mass times the lattice scale, was 0.54(3) in Ref. [|^, as opposed 
to 0.59(6) in Ref. P], giving the ratio of scales, a(/3 = 5.9)/a(/3 = 6.0) = 1.09(13). Quenched 
asymptotic scaling would predict 1.12. 

Let us make several additional points about the data sets shown in Fig.9. First, since 
the physical momentum in these two studies are so different, it is not required that the 
results be in agreement with one another at a given ma value. Second, since the earlier 
study did not use correlated values to monitor the quality of their hts and because of the 
different lattice sizes and (3 values, a comparison between the error bars can not be made 
directly. Ref. Q] does not give chiral charge radii, but we will do a rough extrapolation of 
their implied mM/mp values in the next section to compare with ours. Third, the correlated 
£t shown in Fig. 9 is quite different from what one would expect from an uncorrelated fit of 
the displayed error bars. Note that this is not a special result on the selected correlation 
time interval {t = 14 to 16), but occurs generally in our data whenever the same 3-time-shce 
interval is selected across all k values. 

Experimentally, a monopole form for the spacelike pion form factor gives an excellent £t 
to the data with, however, a pole mass about 4 ± 1% low compared to the mass of the rho 


17|. This value was arrived at from the uncertainty in the £t pole mass given in the hrst 


paper in Ref. |^: 736 ±9 MeV. (The value of the mass in the 1994 Particle Data Group 
review 


ISj is given as 768.1±1.3 MeV.) Note the agreement in the pole mass value in the two 


studies in Ref. O] despite the difference in the studied \q^\ ranges. Thus, it is encouraging 
that the chiral extrapolation of our results agrees in direction and approximate magnitude 
with the small naive vector dominance violation seen in experiments. Of course, this may 
simply be fortuitous since our error bars are about twice the size of the experimental pole 
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mass shift. 


In order to compare with the pion case, we again examine a monopole mass ht for the rho 
meson lowest momentum. Fig. 10 shows the results and Table V lists the mM/fnp ratios for 
the three k values of this study. The linear chiral extrapolation of the data gives 0.73(10). 
The x^d correlated £t is 0.13. It is clear that the error bars here are signihcantly 

larger than for Ref. P]. In addition, although both studies are low compared to vector 
dominance, our hnal extrapolation is about 2.Ter away from the expected vector dominance 
result of mM/mp = 1. Again, note our point above regarding the physical |g^| values of the 
two studies. 

The last quantity we examine is the rho meson magnetic form factor, values of which 
are presented in Table IV. We note that although the conserved lattice current density 
is extended in space, the spatial separations used in the Fourier transform (see Eq.(|TT]) 
above) are unambiguous since the Fourier transform is taken in directions perpendicular to 
the Lorentz index. Unlike the electric form factor, we have no absolute normalization for 
the = 0 value of G'm(g^)- To help evaluate the signihcance of the data, we make the 
assumption that the electric and magnetic form factors have the same dependence for 
our range of momentum transfer. That is, independent of the unknown functional form, we 
simply assume {gp is a constant in g^) 


GM) ®'” 


(30) 


at each value of k, and look at the dependence of Qp on the quark mass. In the nonrelativistic 
limit we expect that gp ~ 2. The nonrelativistic quark model determines the magnetic 
moment in terms of the constituent quark mass, rUg, as 


l^p — 


e 

2mq 


Experimentally, qip = gpel2mp, so that gp = mp/mq, which is rendered less than 2 by the 
binding energy of the rho. See also Ref. for a prediction from a light-front model, in 
excess of 2. 
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In Fig. 11 and Table V we examine the g-factor that results. This quantity seems to 


approach it’s hnal value from below. After doing a linear extrapolation, we hnd that Qp = 
2.25(34). The x\ fhis correlated £t is 0.03. 

In order to estimate the size of the hnite lattice spacing systematic errors in this simu¬ 
lation, we used the lattice spin 0 dispersion relation ^ 


Sink (—) = sink (—) + sin (—), 


and the normalized lattice propgator for a particle of mass tum 


(31) 




sinh\^) -Y.iSin\^^) 
sinK^ {121ft) 


-1 


(32) 


(the analog of Eq. (p^ ) to recalculate the fits in the correlation functions Eqs. ®-(0) above. 
(The dimensionless mass parameter, M, in Eq.(3.17) of Ref. is related to the mass 
defined by the exponential falloff of the two-point function, m by m = 2sinh~^{^■).) Also, 
we made the continuum to lattice replacements [pT|, 


Pi 


E 


Pitt, 


2a ^sin {—), 


2a ^sinh{ 


2 

Ea, 


(33) 

(34) 


for the kinematic factors in these equations. The effect of these replacements was to decrease 
the chirally extrapolated monopole to rho mass ratio, by approximately 3.6% for 

the pion and 2.7% for the rho meson. The extrapolated rho meson g-factor was increased 
by 6.6% by these modihcations. 

To partially investigate the effect of inserting lower momentum, we evaluated the form 
factors at K = 0.152 on lattices of size 20^ x 30. These lattices were also generated at /3 = 6.0. 
They were thermalized by 5000 sweeps and separated by 1000 pseudo-heatbath sweeps. The 
measurements on these larger lattices were done as similar as possible to our 16^ x 24 lattice 
measurements to provide a direct comparison. Using the same scale as before, the lowest 
momentum on these lattices is \q^\ ~ 0.30 GeV^, still about twice that in Ref. [^. The 
dimensionless pion and rho masses measured using 50 conhgurations on f = 15 to 18 hts 
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(t = 0 again labels the time origin) are given by 0.479(4) and 0.547(6), respectively. The 
form factors themselves were again measured with t = 14 to 16 hts on 10 of these lattices. 
We found F = 0.80(3) and Gc = 0.65(8), which gives F/F^^ = 1.04(4) for the pion and 
GcjG^^ = 0.86(11) for the rho meson. When converted into an implied monopole mass, 
these results give = 1.09(11) for the pion and = 0.76(13) for the rho. The 

error bars on the rriM/nip ratios increased signihcantly on the 20^ x 30 lattices both because 
a smaller number of conhgurations were used and because standard error propagation shows 
that the relative error in the monopole mass is SruM/mM = 6F/[2F{1 — F)] in terms of the 
form factor, F. The {1 — F) quantity in the denominator increases the error bars at lower 
momenta. However, comparing with the k = 0.152 results in Tables III and V, we see that 
both rriM/mp ratios have increased. This is what one expects if vector dominance is to be 
reinstated at lower |g^|. 


IV. SUMMARY AND COMMENTS 

We have calculated the pion and rho meson form factors on a 16^ x 24 lattice using the 
methods of current overlap. We limited our measurements to the lowest lattice momentum 
because of a possible violation of continuum dispersion seen in the pion local {E — m) mea¬ 
surement. We have extracted a very clean pion signal which seems to be single exponential 
even near the time edges. In contrast, the rho meson electric signals are relatively noisy and 
seem to be more strongly affected by time-edge effects. We saw that the combined Q(±; t) 
and Q{0;qly,t) values had a strong overlap with the Q{±;qly,t) results, consistent with 
a zero quadrupole form factor. However, because the correct exponential behavior in the 
latter quantity at k = 0.154 and 0.152 was not demonstrated in the data, a signal for the 
quadrupole form factor could not be isolated. We also saw that the rho magnetic correla¬ 
tion functions, Ky{+-, q^, t), develop faster in time than either the pion or rho electric ones, 
reminiscent of a result in Ref. |p . Our error bars on the magnetic correlation function were 
encouraging; assuming Eq.(PD[), we extracted a rho meson g-factor and found g = 2.25(34). 
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This value is consistent with the nonrelativistic quark model value of 2 and with the recent 


results of a light-front model [19 


Three sources of systematic error in our results were investigated. First, by comparing 
pion source time separations of AT = 17 and 23, we saw no statistically signihcant difference 
in the time signals for the lowest momentum. We also found that the lattice was long enough 
to hlter out the correct local {E — m) for both the pion and rho meson. Second, by using 
the lattice spin 0 dispersion relation and dehnitions of kinematical factors, we saw that 
the physical mu/mp ratios for the pion and rho meson were decreased by 3.6% and 2.7%, 
respectively; the rho meson g-factor was increased by 6.6%. In addition, our results on the 
20^ X 30 lattices at k = 0.152 gave larger values of iriM/fTtp than the smaller lattices, but 
with larger error bars. Systematics associated with quenching and scaling have not been 
investigated and must wait future studies. 

A comparison with the pion and rho results in Ref. has been carried out. A simple 
uncorrelated linear extrapolation of the Ref. 0 data displayed in Figs.9 and 10 yields 
a mM/rrip ratio of 1.08(8) for the pion and 0.96(7) for the rho. (The CURFIT routine 
from Ref. was used for the £t and error bars.) Our results for these same quantities 
are 0.97(8) and 0.73(10). If we combine errors in an uncorrelated fashion. Ref. [^ finds a 
rho meson to pion charge radius of 1.13(12), whereas our results give 1.33(21). The results 
in Ref. [§] were extracted at |g^| ~ 0.16 GeV^, whereas here |g^| ~ 0.47 GeV^. Figs.9 
and 10 would seem to indicate a systematic difference in the results at finite ma. 

However, both studies have chirally extrapolated results for the pion which are consistent 
with vector dominance and with one another, within large errors. On the other hand, our 
extrapolated rho meson mu/mp is approximately 2.7a low compared to vector dominance, 
whereas Ref. PJ is consistent with the vector dominance prediction. We noted above that 
our 20^ X 30 lattice results hint that the larger |g^| of the present study may be responsible 
for the differences seen. 

It may help to put our results in the context of measurements of similar quantities in 
the baryon case. We are concentrating on the dimensionless quantity rriM/mp in the present 
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work, similar to the dipole mass ratio mD/rriN that was studied for the nucleon in Ref. 0. 
On the same size of lattice as here, it was found in Ref. [Q that the mu/rriN measurement 
was about 7% low compared to experiment after chiral extrapolation. The extrapolated 
g-factors there were about 10 — 20% smaller in magnitude than experiment. It is reasonable 
to expect there will be similar sized overall systematic errors on the m^/mp ratios and the 
rho meson g-factor in our case. 

Comparison with the experimental results for the pion form factor 0 are suggestive of 
the hnding of a small violation of naive vector dominance in the spacelike pion form factor. 
Experimentally, the the rriM/^p ratio is 0.96(1). Clearly, the error bars on lattice studies 
must be signihcantly reduced before the results can begin to have phenomenological impact 
in this sector. 

The simplest explanation for a small violation of naive vector dominance in the pion form 
factor involves radially excited states of the rho®. Such physics is contained in the lattice 
simulation, even in the quenched approximation. Thus, it is conceivable that the nonzero 
ma results are indicative of this violation at our \q^\. In any case, there must be additional 
physics, outside of naive vector dominance influencing the pion and the rho, resulting in 
their different charge radii. Making stronger contact with these experimental issues will be 
a challenge for lattice QCD in the future and will require signihcantly larger statistics to 
probe closer to the chiral limit and to substantially reduce statistical error bars. 
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Footnotes 


1. New techniques such as in S. J. Dong and K. F. Liu, Phys. Lett. B328 (1994) 130, 
and Y. Kuramashi et. ah, Phys. Rev. Lett. 72 (1994) 3448, are beginning to change 
the situation for the better. 


2. We are grateful to T. Draper for telling us about this technique. 

3. Note that it is always possible to £x a single spatial point in any n-point correlation 
function without introducing momentum smearing. 


4. 

5. 


Such an assumption is implicit also in the results of Ref. @ for the rho meson. 


See for example Ref. where two additional 
1420 and 1770 MeV are sufficient to explain the 
energy region 1.35 < ^/s < 2.4 GeV. 


vector mesons with masses of 
timelike pion form factor data 


about 
in the 
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TABLES 


TABLE 1. Hadron masses, xl gives the chi-squared per degree of freedom for the fit. 


Particle 

K = 0.154 

0.152 

0.148 

pion 

0.366(10) 

0.479(8) 

0.676(6) 

xl 

.45 

.43 

.28 

rho 

0.46(2) 

0.55(1) 

0.718(8) 

Xd 

.23 

.19 

.05 

TABLE II. Pion form factors for the lowest spatial momentum. 

indicates the vector 

dominance value for the lattice value of 

nip from Table I. Xd gives the chi-squared per degree of 

freedom for the fit. 




Quantity 

F = 0.154 

K = 0.152 

F = 0.148 

F 

0.58(4) 

0.65(2) 

0.74(3) 

F 

-pVT7 

0.91(7) 

0.94(4) 

0.94(3) 

xl 

0.65 

0.20 

1.02 

TABLE III. Values of monopole mass divided by rho meson mass, mM/'^ip for the lowest 

momentum fit to the 

pion form factor, k 

c indicates the correlated chiral extrapolation. 

Kc 

K = 0.154 

K = 0.152 

F = 0.148 

0.97(8) 

0.89(6) 

0.91(5) 

0.89(6) 
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TABLE IV. Gc and Gm are the charge and magnetic form factors of the rho meson. G^^ indi¬ 
cates the vector dominance value for the lattice value of rup from Table I. Xd gives the chi-squared 
per degree of freedom for the fit. 


Quantity 

K = 0.154 


K = 0.152 

K = 0.148 

Gc 

0.46(6) 


0.55(5) 

0.68(4) 

Gc 

Gvu 

0.74(9) 


0.81(7) 

0.87(6) 

xl 

0.43 


0.33 

0.23 

Gm 

0.93(15) 


1.07(9) 

1.14(7) 

xl 

0.79 


0.83 

0.79 

TABLE V. ruM/' 

nip is the ratio of monopole mass 

to rho meson mass 

obtained from the Gc 

value for the rho meson. The g-factor 

assumes Eq.(30). 



Quantity 

Kc 

K = 0.154 

K = 0.152 

K = 0.148 

niM/mp 

0.73(10) 

0.73(8) 

0.75(7) 

0.77(8) 

g-factor 

2.25(34) 

2.03(38) 

1.93(21) 

1.67(10) 
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Figure Captions 


1. (a) Graphical representation of the four-point function being calculated. The time 
positions of the particle interpolation helds (shown as shaded circles) are fixed, but 
the charge and current densities, located on different flavor lines and symbolized by 
the “X”, may be time separated, (b) Example of a disconnected (sea quark) contribu¬ 
tion. The disconnected quark loop is actually connected to the quarks by an arbitrary 
number of gluon lines. 

2. (a) Fourier transform of pion charge overlap function at zero relative time separation, 
Q{q‘^,0), for smeared-to-smeared interpolation helds (k = 0.154) as the charge den¬ 
sities are moved together in time between the sources, (t = 0 labels the hrst time 
slice of the lattice.) (b) Fourier transform of the same quantity for smeared-to-point 
interpolation helds. (The point interpolation held, located at t = 24, is unsummed.) 

3. Comparison of pion overlap correlation function, Q{q‘^,t), versus relative time separa¬ 
tion, t, for a time separation AT = 17 (squares) between the interpolation helds and 
a time separation of AT = 23 (diamonds). Both calculations use smeared-to-smeared 
quark propagators with k, = 0.154. 

4. Local {E — m) measurements for the pion overlap function, Q{q‘^,t), versus relative 
time separation, t, compared with continuum dispersion for (a) \qmin\ = f and (b) 

= V^\qmin\- Squares are k = 0.154, triangles are k = 0.152, and circles are 
K = 0.148 results. Continuum dispersion is indicated by the solid lines; the associated 
spin 0 lattice dispersion relation is indicated by dotted lines underneath. 

5. Correlation functions, Q{q‘^,t), for the lowest spatial momenta (|^ = |) of the pion 

versus relative time separation, t. Results are at k = 0.148 (circles), k = 0.152 (tri¬ 
angles) and K = 0.154 (squares). (The kinematical factor has been removed 

from Q{q^, t).) The best uncorrelated single exponential hts on the 14-16 time interval 
are shown. 
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6. The electric correlation function (5(±; Qz, t) for the rho meson for the lowest spatial 

momenta (|^ = f) versus relative time separation, t. Results are at k = 0.148 (circles), 
K = 0.152 (triangles) and k = 0.154 (squares). (The kinematical factor has 

been removed from The best single exponential hts on the 14-16 time 

interval are shown for k = 0.148 and 0.152; the best £t for the 13-15 interval is shown 
for K = 0.154. 

7. The magnetic correlation function Ky{+] t) for the rho meson for the lowest spatial 

momentum (|^ = |) versus relative time separation, t. Results are at k = 0.148 
(circles), k = 0.152 (triangles) and k = 0.154 (squares). (The factor iq^ has 

been removed from Ky{+',qx,t).) The best single exponential fits on the 7-9 time 
interval are shown. 

8. The quadrupole correlation function Qp{±;qly,t) for the rho meson for the lowest 

spatial momentum (|^ = |) versus relative time separation, t. (The kinematical 
factor has been removed.) Results are at k = 0.148 (circles), k = 0.152 

(triangles) and k = 0.154 (squares). No fits of these functions were made. 

9. Monopole mass divided by the rho meson mass, rriM/'^p, for the pion as a function of 
dimensionless quark mass (ma) as determined from the lowest momentum measure¬ 
ments (squares). The darkened square gives the chiral extrapolation. The triangles 
give results extracted from Table 1 of Ref. P] . 

10. Monopole mass divided by the rho meson mass, itim/^pi for the rho meson as a 
function of dimensionless quark mass (ma) as determined from the lowest momentum 
measurements (squares) assuming the quadrupole form factor is zero. The darkened 
square gives the chiral extrapolation. Triangles are the results extracted from Table I 
of Ref. §. 

11. g-factor for the rho as a function of dimensionless quark mass (ma) assuming Eq.(^). 
Darkened square gives the chiral extrapolation. 
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